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Abstract 

Using BFS, a new semiempirical method for alloys, we study the surface structure of 
fee ordered binary alloys in the LI 2 structure (NiaAl and CU 3 AU). We show that the 
surface energy is lowest for the mixed- composition truncation of the low-index faces of 
such systems. Also, we present results for the interlayer relaxations for planes close 
to the surface, revealing different relaxations for atoms of different species producing a 
rippled surface layer. 
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1. Introduction 


With the recent widespread use of semiempirical methods, as well as the ever-increasing 
computational capability, there has been an enormous increase in the research in computa- 
tional material science that, until not long ago, was considered intractable. Surface structure 
determinations also saw a tremendous development with the advent of new experimental 
techniques as well as new theoretical approaches for elemental crystals as well as alloys. 
Recently, much effort has been devoted to the study of surface relaxation in metals and, to 
a lesser extent, alloys. Several recent experiments have provided insights in the phenomena 
of surface relaxation and composition, in the case of alloys, and correspondingly a number 
of theoretical studies have shown good general agreement with experimental results. 

However, there is still a great deal of uncertainty in certain areas, due to limitations 
inherent in experimental techniques and also to the lack of alternative studies to verify 
previous results. Semiempirical methods are particularly useful tools, providing an efficient 
and economical way of investigating problems and the consistency of conclusions drawn 
from experiment. 

The wealth of experimental studies of surface relaxation on pure metallic surfaces is 
not matched for alloys [1]. However, in spite of the small number of experimental studies 
[2,3,4], there seems to be slow but sure progress in the field, as the available theoretical 
tools for modelling become more accurate. 

The first experiment, in 1984, that provided detailed information on the atomic po- 
sitions of surface atoms in a truncated ordered alloy is the low-energy electron diffraction 
(LEED) intensity analysis of Davis and Noonan of a NiAl(llO) surface [2]. They found 
strong evidence for a rippled surface, where the A1 sites of the top layer (in the mixed- 
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composition truncation) are displaced above the Ni sites by approximately 0.22 A. This 
result was quickly followed by the calculation performed by Chen, Voter and Srolovitz us- 
ing the embedded atom method (EAM) [5], which confirmed the main features found in 
the experiment. EAM was later used to investigate similar phenomena in other ordered 
alloys: Foiles and Daw presented a complete study of Ni 3 Al (Ll 2 structure) [6], followed by 
Foiles work on ordered surface phases of Au on Cu [7], and Lundberg’s extensive study of 
surface segregation and relaxation of Pt-Ni alloys [8]. At the same time, new experimental 
LEED results on Ni 3 Al were reported by Sondericker and coworkers [3], finding a similar 
rippled structure in Ni 3 Al (100) faces. Finally, a low-energy ion scattering spectroscopy 
(LEISS) experiment by Wang and coworkers provided similar data for the Cu 3 Au system, 
a classic ordering alloy [4]. Their work followed the LEISS results concerning the surface 
composition of the top atomic layers [9]. This system was also the subject of a very recent 
study by Wallace and Ackland using a molecular statics algorithm with Finnis- Sinclair (FS) 
many-body potentials [10]. 

The purpose of this work is to provide new results concerning the above mentioned sys- 
tems, using a new semiempirical method for alloys recently developed by Bozzolo, Ferrante 
and Smith (BFS) [11], which has shown great promise in previous applications to problems 
of alloy structure [12]. Due to the simplicity of the method and the ensuing computational 
efficiency, this application of BFS to the surface structure of ordered alloys can be easily 
extended to a number of systems for which there are no other theoretical or experimental 
studies to date [12]. 

2. The BFS method 

The BFS method is based on the idea that the energy of formation of an arbitrary 
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alloy structure is the superposition of individual contributions e, of non- equivalent atoms 
in the alloy [12-13]: 

e i = e f + 9i( € ? ~ e f°)- (1) 

e, has two components: a strain energy e 5 that accounts for the actual geometrical distri- 
bution of the atoms surrounding atom i , computed as if all its neighbors were of the same 
atomic species, and a chemicals nergy e^, which takes into account the fact that some of the 
neighbors of atom i may be of a different chemical species. The coupling function g t , ensures 
the correct asymptotic behavior of the chemical energy contribution. The strain energy of 
a pure defect crystal is calculated assuming that every neighbor of atom i is of the same 
species X. ef is then computed with any technique (first-principles methods, semiempirical 
techniques, etc.). For ef we interpret the chemical composition as a defect of an otherwise 
pure crystal. We represent this defect by ‘perturbing’ the electronic density in the overlap 
region between dissimilar atoms and locating them at equilibrium lattice sites of atom i. 
The ideas of equivalent crystal theory [13] are used to develop a procedure for the evaluation 
of the energy associated with this ‘defect’. To free the chemical energy of structural defect 
energy which should only be included in the strain energy, we reference ef to a similar 
contribution where no such perturbation is included (ef°). Finally, the coupling function 
< 7 , is defined as gi = e~ a ' , where af is a solution of ef = — E l c |l — (1 4- af )exp{— af )j , 
and where E' c is the cohesive energy for atom i. We direct the reader to ref. 12 where 
a detailed description of the calculation of the strain and chemical energy contributions is 
provided. Except for two parameters determined by fitting to experimental or theoretical 
alloy properties, the method relies on pure element properties. Moreover, the pairwise char- 
acter of the interaction between dissimilar species facilitates application to multicomponent 
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systems with no further experimental or theoretical input. Within the framework of BFS, 
the calculation of defect energies requires information on the atomic positions only. The 
input parameters (see ref. 12 for a complete list) are readily available for a variety of fee 
and bcc alloys. 

3. Multilayer relaxation of pure crystals 

Before proceeding to the calculation of multilayer relaxation in alloys, we will discuss 
some features of theoretical calculations of these quantities. Ref. 1 provides a reasonably 
large sample of both experimental and theoretical results for changes in interlayer spacing 
in pure crystals. In all cases, the theoretical techniques used rely either on the use of input 
data (generally experimentally determined) or on certain approximations for some of the 
variables of relevance. Necessarily, the results will depend on such choices. Multilayer 
relaxations involve at best very small changes in position, and correspondingly, comparable 
changes in surface energy, whose minimization is the criterion used to determine the final 
interlayer spacings. Thus, the search for a minimum of the surface energy, as accurate as the 
minimization technique might be, will be strongly influenced by the two factors indicated 
above: the approximations used and the shallowness of the minimum in the surface energy 
surface resulting from small changes in the input parameters. As a consequence, to quote 
just one value for each of the changes in interlayer spacings as is ordinarily done, might 
not reflect the ambiguities in these calculations. In this paper we adopt a different path: 
to each theoretical prediction, we will attach an estimate of the possible errors due to 
any of the reasons mentioned above. Although there is no certain way to determine such 
errors (after all, the predictions are, within their own framework, exact), we will see that 
changes on the order of 1 % in the surface energy can generate quite interesting variations 
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in the relaxation schemes predicted. In particular, within the framework of ECT, such 
small changes in the surface energy can be easily obtained by changing any of the input 
parameters (lattice constant, cohesive energy, bulk modulus) by a similar amount, well 
below the usual experimental errors in the determination of such quantities. 

To illustrate this issue, we will focus our attention on the surface structure of some 
fee pure metals (Al, Au, Cu and Ni). As can be seen in Tables 2-11 of ref. 1, previous 
theoretical and experimental studies show a wide spread in the predictions of the changes 
in interlayer spacings for the (100) and (110) surfaces. Even results obtained within the 
same theoretical technique (EAM, ECT) do not agree with each other (due to different 
fitting procedures of the embedding function in the case of EAM and different input data 
in both cases). Although there is general qualitative agreement, regarding the contraction 
or expansion pattern found for successive layers, in some cases the theoretical values show 
poor agreement with experimental results (see, for example, Al (100)). The ECT results 
(from refs. 1 and 13) also highlight this inconsistency. The difference between the values 
obtained in this work and those from previous applications of ECT is easily traceable to 
slightly different values of some of the input parameters. 

4. Multilayer relaxation of ordered alloys. 

As mentioned above, in order to account for these and other ambiguities in the calcu- 
lation, we investigated the change in predicted relaxations due to small changes in the rigid 
surface energy. We thus defined ‘error bars’ in such way that all the intermediate values 
so obtained predict variations in surface energies within that tolerance. Needless to say, 
this range of values does not include all the possible sets (Adi 2 , Ac?23) that correspond to 
surface energies within the allowed values. It is interesting to note, however, that in most 
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cases, all the experimental as well as theoretical predictions fall within the range of values 
obtained in this fashion. 

It should be noted that when comparing our theoretical predictions with available 
experimental results, the error bars quoted in each case are not rigourously comparable. 
However, we choose to do so with the only purpose of giving a complete description of the 
results obtainable with EOT (for pure crystals) and BFS (for alloys), once uncertainties 
in the input parameters are taken into account. To illustrate this point, we first discuss 
the surface energies and multilayer relaxations of the unreconstructed low-index surfaces of 
pure Al, Ni, Cu and Au crystals. In Table 1 we display the ECT predictions for the surface 
energies and compare the results with typical experimental values for polycrystalline samples 
[14,15]. The agreement is excellent in all cases. We note that experimental values for the 
surface energies are for polycrystalline surfaces, thus could be strongly dominated by the 
predominant surface plane. 

In table 2 we compare results for the multilayer relaxations of the first two interlayer 
spacings for those cases for which recent experimental data is available [16-24]. Once again 
the agreement is excellent, as it was shown in previous ECT studies of surface structure 
[1], The inclusion of the theoretical ‘error bar’, as mentioned above, allows for a better 
comparison with experiment as it shows that for most cases, small changes in the input 
parameters of the method suffice to account for the whole range of possible experimental 
results. The exceptions are Al(100) and Al(lll), where the outward relaxation of the surface 
layer has been attributed to an electron promotion effect [25]. Semiempirical methods 
(ECT, EAM. etc.), unless specifically designed to do so, do not generally allow for such 
fine electronic structure effects, thus it is not surprising that our results for Ad 12 in these 
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cases predict surface layer contractions, even when the ‘error bar’ is taken into account. 
For completeness we also include results for the surface relaxation when only that plane is 
allowed to relax, in order to single out the influence of subsequent interlayer spacing changes 
on the surface plane. Again, the agreement with available experimental data is very good 
in all cases. 

As mentioned above, there are few theoretical or experimental studies of ordered alloy 
surfaces [2-10]. In this paper, we focus our attention on two cases (Ni 3 Al and Cu 3 Au, in the 
LI 2 structure) which have been the subject of recent studies [3-6]. A complete presentation 
of the corresponding results for a larger number of ordered structures as well as different 
binary alloys of both bcc and fee elements, for which no theoretical or experimental data 
exists, will be published elsewhere. 

a. Cu 3 Au: Table 3 displays the results for the unrelaxed and relaxed surface energies 
(in ergs/cm 2 ) as obtained with BFS and with FS many-body potentials [10]. Both methods 
predict, as expected, lower surface energies for the mixed-composition (100) and (110) 
truncations. This feature has been experimentally proven via a low-energy ion scattering 
study which detected equal parts of Cu and Au in the top layer [9]. ECT and FS results 
also agree on the relative change in surface energy once the top-most layers are allowed to 
relax, in spite of the fact that the FS values are 50 % smaller than the BFS ones. As is also 
to be expected, the surface energies of (100)1:0 and (110)1:0 faces are comparable to the 
corresponding values for single Cu crystals. The corresponding relaxations are quoted in 
table 4. In order to avoid ambiguities in determining the exact atomic positions from the 
entries in table 4, we present the relaxations as the percentage change in interlayer spacing 
from the unrelaxed case to the one measured from the relaxed position to the unrelaxed 
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location of the plane immediately below. We also include the BFS predictions for the pure 
Cu truncations of the ( 100 ) and ( 110 ) planes. Although it is to be expected that the top 
layer relaxation will change as deeper layers are allowed to relax, any ensuing changes would 
be small, not affecting the conclusions drawn from our results. 

For the Cu 3 Au ( 100 ) 1:1 Cu:Au case, the results in table 4 imply a rippling of 0 . 148 ± 0.023 
A, which amounts to 3 . 97 ± 0.62 % of the lattice parameter determined for this alloy ( 3.73 
A). This result compares very well with the 3.77 % rippling (A1 out, Ni in) obtained using 
FS potentials [ 10 ]. A similar situation is found for the ( 110 ) 1:1 Cu:Au surface, where we 
find the rippling to be 0.135 ± 0.04 A ( 3 . 6 ± l.l % of the lattice parameter), whereas FS 
potentials predict a rather smaller change of 1.9 %. For the mixed-composition ( 111 ) 3:1 
Cu:Au surface, BFS predicts a rippling of 4 . 6 ± 0.4 % thus agreeing with FS results and 
experimental evidence that the Au atoms are farther out than the neighboring Cu atoms 
in mixed-composition surfaces. 

b. Ni 3 Al. The surface energies of relaxed ( 100 ), ( 110 ) and ( 111 ) surfaces are shown 
in Table 5 , where we compare our results with the EAM study of Foiles and Daw [ 6 ]. As 
found for the Cu 3 Au case, the mixed-composition truncations always have a lower surface 
energy. The differences between the EAM and BFS predictions are consistent with previous 
calculations for pure metals, where the EAM results are typically 50 % lower than the 
experimental ones. Surface relaxations are indicated in Table 6, using the same format and 
notation of table 4. From these results we extract the following values for the gap between 
Ni and A1 atoms in the mixed-composition ( 100 ), ( 110 ) and ( 111 ) surfaces: 0 . 12 ± 0.04 A, 
0 . 07 ± 0.05 A and 0 . 16 ± 0.03 A, respectively. A similar trend, but with somewhat smaller 
values for the rippling are obtained from EAM [ 5 ]: 0.09 A, 0.06 A and 0.07 A. A different 
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EAM calculation [7] predicts a 0.06 A separation between Ni and A1 atoms in all three 
surfaces. Recent LEED data [3] show Adj 2 (Ai) = —2.73% and Ac£i 2 (Al) ~ 0 (i.e., A1 is 
displaced outward with respect to Ni) and a rippling of 0.02±0.03 A for the (100) surface 
In conclusion, we have applied BFS for the study of multilayer relaxation of specific 
ordered alloy surfaces for those cases where comparison with experiment and other theo- 
retical approaches is possible. We found our Cu 3 Au results in agreement with experiment 
and other calculations in that the mixed layer termination had a lower surface energy than 
the Cu terminated layer and that surface relaxation studies indicated that there would be 
rippled surfaces with a preferred outward relaxation for the Au atoms. For Ni 3 Al similar 
agreement was obtained for the surface energy as a function of surface composition as well 
as for the outward relaxation of A1 atoms on the surface. Finally, we have pointed out 
that considerable caution must be exercised in presenting quantitative results for surface 
relaxations in that small changes in input parameters can cause substantial quantitative 
changes due to shallowness of the energy surface. 
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Technique 

Al 

Cu 

Ni 

Au 

Exp. [14] 

1200 

1790 

2270 

1560 

Exp. [151 

1140 

1780 

2380 

1500 

Exp. [15] 

1180 

1770 

2240 

1540 

ECT(IOO) 

1203 

2309 

2982 

1546 

ECT(llO) 

1284 

2373 

3073 

1621 

ECT(lll) 

856 

1767 

2274 

1136 


Table 1: Surface energies (in ergs/cm 2 ) 


Element 

Face 

Experiment 


ECT 

ECT 

(two-layers) 



Adj2 

Ac?23 

Ref. 

Ad\2 

Adi2 



(100) 

+ 1.8 


[16] 

-4.68+1.62 

-5.05+1.58 

+3.35+0.80 

Al 

(110) 

-8.5+1.0 

+5.5+1. 1 

[17] 

-8.29+2.35 

-9.53+3.58 

+ 1.90+2.24 


(111) 

+ 1.7+0.3 

+0.5+0. 7 

[18] 

-3.67+1.21 

-3.94+1.19 

+2.75+0.61 


(100) 

-3.2+0.5 


[19] 

-3.53+1.68 

-3.82+1.68 

+2.48+0.85 

Ni 

(110) 

-9. 0+1.0 

+3.5+1. 5 

[20] 

-6.32+2.44 

-6.55+3.63 

+0.34+2.24 


(111) 

-1.2+ 1.2 


[21] 

-2.89+1.29 

-3.10+1.25 

+2.12+0.63 


(100) 

-2.1 

+0.45 

[22] 

-3.52+1.74 

-3.81+1.70 

+2.47+0.86 

Cu 

(110) 

-7. 5+1. 5 

+2. 5+1. 5 

[23] 

-6.31+2.46 

-6.51+3.83 

+0.29+2.44 


(111) 

-0.7+0. 5 


[24] 

-2.88+1.30 

-3.10+1.25 

+2.12+0.63 


Table 2: Surface relaxations of Al, Cu and Ni as percentages of the rigid interplanar spacings. 


Face 

Finnis-Sinclair 

BFS 


U nrelaxed 

Relaxed 

Unrelaxed 

Relaxed 

(100)1:1 

896 

865 

2119 

1810 

(100)1:0 

1192 

1171 

2478 

2247 

(110)1:1 

1051 

1024 

2397 

2337 

(110)1:0 

1240 

1173 

2873 

2699 

(111)3:1 

882 

863 

1626 

1577 


Table 3: Unrelaxed and relaxed surface energies (in ergs/cm 2 ) of Cu 3 Au. 
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Layer 

Atom 

(100)1:1 

(100)1:0 

(110)1:1 

(110)1:0 

(111)3:1 

1 

Cu 

-2.12±0.41 


-2.45±1.14 

-3.70±1.02 

+ 14.35±0.26 


Au 

+5.82±0.81 

- 


- 

+22.29T0.53 

2 

Cu 

+5.61±0.37 

+ 12.87±0.50 

+10.14±1.37 

+12.14±1.37 

+21.20±0.24 


Au 

- 

+3.23±1.23 

- 

-0.70±1.30 

+2.23±1.92 


Table 4: Surface relaxations of Cu 3 Au Ll 2 surfaces. See text for definition of the percentage 
change. 


Face 

EAM [6] 

BFS 

(100)1:1 

1620 

2852 

(100)1:0 

1885 

3168 

(110)1:1 

1730 

3117 

(110)1:0 

1920 

3964 

(111)3:1 

1645 

2411 


Table 5: Surface energies of Ni 3 Al (in ergs/cm 2 ). 


Layer 

Atom 

(100)1:1 

(100)1:0 

(110)1:1 

(110)1:0 

(111)3:1 



-1.33±0.70 

+3.41±0.72 

-4.63±1.65 

-5.03±1.31 

+8.56±0.55 



+5.25±1.60 

- 

+0.80±2.85 


+ 16.28±1.40 

2 

Ni 

+4.59T0.63 

+8.82±0.81 

+3.55±2.19 

+5.19±2.40 

+ 13.69±0.53 


A1 

- 

+3.39±1.46 

- 

+0.92±2.74 

+4.04±1.53 


Table 6: Surface relaxations of Ni 3 Al Ll 2 surfaces. See text for definition of the percentage 
change. 


14 















REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

Public reporting burden tor this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate tor Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington. VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 

1. AGENCY USE ONLY {Leave blank) 

2. REPORT DATE 

December 1993 

3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 

4. TITLE AND SUBTITLE 

Multilayer Relaxation and Surface Structure of Ordered Alloys 

5. FUNDING NUMBERS 

WU-505-90-53 

6. AUTHOR(S) 

Robert J. Kobistek, Guillermo Bozzolo, John Ferrante, 
and Herbert Schlosser 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

E-8275 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESSES) 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


NASA TM-1 06433 


11. SUPPLEMENTARY NOTES _ 

Robert J. Kobistek, Victoreen, Inc. 6000 Cochran Road, Solon, Ohio 44139-3395; Guillermo Bozzolo, Analex Corporation, 3001 
Aerospace Parkway, Brook Park, Ohio 44142, (work funded by NASA Contract 3-25776); John Ferrante, NASA Lewis Research 
Center; and Herbert Schlosser, Cleveland State University, Physics Department, Cleveland, Ohio 441 15. Responsible person, John 
Ferrante, (216) 433-6069. 

12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category 26 


12b. DISTRIBUTION CODE 


13. ABSTRACT (Maximum 200 words) 

Using BFS, a new semiempirical method for alloys, we study the surface structure of fee ordered binary alloys in the 
Ll 2 structure (Ni 3 Al and Cu 3 Au). We show that the surface energy is lowest for the mixed-composition truncation of 
the low-index faces of such systems. Also, we present results for the interlayer relaxations for planes close to the 
surface, revealing different relaxations for atoms of different species producing a rippled surface layer. 


14. SUBJECT TERMS 

Alloys; Surface gold; Copper; Nickel; Aluminum 


15. NUMBER OF PAGES 

16 


16. PRICE CODE 

A03 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Sid. Z39-18 
298-102 
















